Molecular Plant-Microbe Interactions®

Accessory genes are thought to contribute to fungal adaptation and pathogenicity by modulating host immunity, while core genes play crucial roles in maintaining fundamental biological processes. Rust fungi (order Pucciniales) are obligate biotrophic plant-pathogens and infect economically relevant crops. Here, we characterize core and accessory gene repertoires across rust fungi, with a particular focus on Phakopsora pachyrhizi, the causal agent of Asian soybean rust. Across Pucciniales genomes, accessory genes represented the largest fraction of gene content (~44.6% on average), whereas core genes accounted for a smaller proportion (~18-35%). Notably, variations in accessory gene content among rust fungi are perhaps attributed to lineage-specific gene expansions and losses. Core gene content was positively correlated with total gene number across Pucciniales genomes, suggesting retention after gene duplication events, consistent with their essential biological functions. Among P. pachyrhizi genes expressed during soybean infection, core effectors were associated with cysteine-rich proteins, pectin-degrading enzymes, and SPFH/Band 7 family, while accessory effectors included phosphatidylethanolamine-binding proteins, trehalose phosphatases, and CFEM domain-containing proteins. The in-plant induced core and accessory genes in P. pachyrhizi also comprised multiple families of CAZymes (GH5/GH7 cellulases, CE5 cutinases, CE8 pectinesterases, CE4/GH18 chitin-modifying enzymes); proteases (aspartyl proteases, serine carboxypeptidases, alpha/beta hydrolases); transporters (amino acid permeases, ferric reductase-like transmembrane proteins, and OPT oligopeptide transporter), and transcription factors (bZIP, GATA zinc finger, STE-like, and homeobox KN). Our study highlights that core and accessory gene families have shaped P. pachyrhizi-soybean interactions, identifying promising targets for functional studies aimed at elucidating host-adaptation mechanisms in rust fungi.

Multiple recognition events of pathogen-secreted effectors by immune receptors confer robust disease resistance in plants. Understanding the underlying mechanisms facilitates the discovery and deployment of valuable resistance genes for crop protection. Xanthomonas euvesicatoria causes devastating bacterial spot disease in solanaceous crops but cannot infect the wild relative Solanum americanum. Here, we identified X. euvesicatoria type III effectors (T3Es) that induce cell death in S. americanum when transiently expressed in leaf. By quantifying immune responses of the S. americanum SP2273 accession to X. euvesicatoria multiple-T3E knockout mutants, we demonstrated that at least nine T3Es (AvrBs2, XopAP, XopAU, XopE1, XopJ3, XopM, XopN, XopX, and XopZ1) collectively contribute to effector-triggered immunity (ETI). Among these, AvrBs2 and XopJ3 were the primary drivers of ETI, eliciting robust cell death and defense gene expression when naturally delivered into plant cells. We next generated S. americanum lines concomitantly edited at the corresponding immune receptor loci, SaBs2 and SaZAR1 (SP2273-bz). Genetic complementation of SP2273-bz confirmed that AvrBs2 and XopJ3 are specifically recognized by each of the four SaBs2 homologs and by SaZAR1, respectively. Moreover, enhanced growth of X. euvesicatoria on the characterized SP2273-bz line indicated that both effectors are required for bacterial multiplication in S. americanum. Together, our findings provide a framework for understanding the mechanisms of ETI-mediated resistance and establish a genetic foundation for resistance breeding in solanaceous crops.

Botrytis cinerea is a necrotrophic plant pathogen with an extremely wide host range. During invasion, the fungus induces rapid host cell death and proliferates in the necrotic tissue. The mechanisms of host killing are still incompletely understood, they involve secretion of lytic enzymes, phytotoxic metabolites and cell death inducing proteins (CDIPs). We have previously shown that the sequential knockout of up to 12 CDIPs leads to a substantial reduction of virulence of B. cinerea mutants. In this study, we have identified additional CDIPs and generated an extended mutant series culminating in a 29x mutant which is deficient in most currently known CDIPs and two phytotoxic metabolites. These mutants are strongly reduced in infection, but still induced necrosis and grey mould symptoms, demonstrating that additional determinants of host killing remain unidentified. Overexpression of the highly phytotoxic Nep1 in a 22-fold CDIP mutant failed to increase its virulence. Reevaluation of several CDIPs previously described as virulence factors revealed for most tested CDIPs no major contribution to pathogenesis. Together with the observation that none of the CDIPs are specific to B. cinerea, our data question a particular role of CDIPs for necrotrophic pathogenesis. In contrast, generation of a mutant lacking all six predicted endo-polygalacturonases confirmed their major but not exclusive role for tissue degradation and infection.

The contribution of accessory or conditionally dispensable chromosomes to host-specific virulence was first demonstrated in members of the Fusarium solani species complex (FSSC) that are pathogens of garden pea, Pisum sativum L. The phenomenon has since been shown to exist in many fungal plant pathogens, including the closely related F. oxysporum species complex (FOSC). Genome analysis of members of the FSSC and FOSC pathogenic on pea revealed a diverse size range of the accessory genome of these fungi. Despite the ~65 million years of diverging time, regions on a chromosome known to carry host-specific virulence factors for pea, including the cytochrome P450 pisatin demethylase (PDA) and other pea pathogenicity (PEP) genes, were present in all genomes of these pea pathogens. Genes directly involved in virulence on pea - PEP2, PDA, and PEP5- were the most frequently clustered together. Transcriptome analysis of fungal mycelia treated with the pea phytoalexin pisatin, identified 1,155 differentially expressed genes where many were involved in cellular stress responses. As wilt pathogens that invade host xylem, members of the FOSC encode more putative effectors, when compared to those in the FSSC, and several FOSC effectors were identified to confer race specificity. The conservation of part of the accessory genomes across two evolutionarily diverged species complexes suggests a common origin. Horizontal transfer of accessory chromosomes containing genetic loci involved in pathogenesis for garden pea offers a parsimonious explanation of the polyphyletic origin of host specificity.

O_LIThe type II secretion (T2S) system is conserved across the Xanthomonas lineage, yet its contributions to pathogenicity and secreted protein repertoires are poorly defined. We demonstrate that T2S systems in Xanthomonas pathovars with divergent hosts and lifestyles are required for disease. C_LIO_LIIn planta quantification of cell wall compositional changes during infection by Xanthomonas euvesicatoria (Xe) revealed that T2S-dependent depletion of galacturonic acid occurs during host colonization, providing experimental evidence for T2 effector (T2E)-mediated cell wall remodeling. C_LIO_LIUsing an in planta label-free proteomics approach, we identified two known and 20 new Xe T2Es from tomato apoplast, many with annotated functions in polysaccharide and protein cleavage. Growth assays on plant cell wall extracts and purified substrates revealed T2S-mediated metabolization of plant cell wall polysaccharides and proteins not only by Xe, but also by Xanthomonas axonopodis pv. glycines (Xag) and Xanthomonas campestris pv. campestris (Xcc). Interestingly, comparative sequence analysis revealed that the T2E repertoires have diversified among these pathogens, with differences in protease repertoire being the most pronounced. C_LIO_LIOur methodology establishes a framework for T2E discovery, enabling future functional dissection of this understudied effector class and its crosstalk with other bacterial virulence factors. C_LI

A highly contiguous and complete reference genome of Cercospora cf. flagellaris, the causal agent of foliar disease on many plant hosts including Cercospora leaf blight of soybean, was assembled using a combination of PacBio and Illumina sequencing reads. The genome assembly is 33.72 Mb in length and consists of 14 nuclear scaffolds and one mitochondrial contig. Four scaffolds have telomeric repeats on both ends and represent fully assembled chromosomes, while nine scaffolds represent partially assembled chromosomes with telomeric repeats on one end. The assembly has an N50 of 2.90 Mb and an L50 of 5 scaffolds. Genome annotation identified 11,268 genes, of which 947 and 360 were predicted to encode secreted proteins and effectors, respectively. Additionally, 512 genes were predicted to encode carbohydrate-active enzymes and 60 biosynthetic gene clusters were annotated. Taken together, this annotated genome assembly will be a valuable resource for genomics, host-pathogen interactions, and population biology research in this economically important pathosystem.

The Pseudomonas syringae species complex includes major crop pathogens that use a type III secretion system (T3SS) to inject effectors into plant cells, suppressing immunity and promoting disease. The cherry canker pathogen Pseudomonas amygdali pv. morsprunorum (Pam) carries the effector gene hopAR1 on a prophage, PamPP1, which belongs to a novel Caudoviricetes family widespread across the P. syringae complex and likely acquired before pathovar divergence. Deletion of PamPP1 shows that this prophage enhances Pam virulence independently of hopAR1, instead it alters the T3SS operon expression both in vitro and in planta. These prophage-driven transcriptional changes likely reshape how Pam interacts with plant immunity, highlighting how bacteriophages rewire bacterial transcriptomes and contribute to the evolution and emergence of plant diseases.

The fungal pathogen Pyricularia oryzae is notorious for causing blast disease in various important cereal crops, including wheat, rice, millet, and oat. Whole-genome-informed data on this pathogen are necessary to better understand the host adaptability of the fungus, including identifying key determinants of infection to enable more precise disease control. Here, we report highly contiguous genome sequences (using long-read PacBio technology) of two isolates from rice paddies in Portugal, M22.7 and T22.2, which exhibit distinctly aggressive symptoms in rice. Both mitochondrial and nuclear sequences were characterised in this study. The resulting nuclear genomes have assembly lengths of 46.4 Mb for M22.7 (198x coverage) and 46.3 Mb for T22.2 (163x coverage), with near-complete BUSCO completeness (98.8%) and a 0% contamination score (EukCC). Phenotypic analysis showed M22.7 to be more virulent than T22.2, which may be explained by the lower number of predicted effector genes and higher transposable element content in M22.7 relative to T22.2. This announcement represents the first genome resource for natural isolates of P. oryzae from Portugal in over 20 years, filling an important data gap from a major European rice-producing country that produces locally adapted rice varieties under specific agro-environmental conditions (near the Atlantic coast).

Potato leafroll virus (PLRV) is an economically important viral disease of potato (S. tuberosum). Genetic resistance to this phloem-limited virus is rare, and no cloned resistance (R) genes have been reported. Rladg confers resistance to PLRV in an Andean potato landrace, LOP-868 (Velasquez et al. 2007). We identified the functional Rladg gene as a homolog of the tomato TIR-NLR-encoding Bs4. Rladg interacts with the serine protease domain of the PLRV protein P1, which is essential for virus replication. This recognition is independent of the proteases enzymatic activity, and the Rladg immune receptor oligomerizes upon direct association with the protease. Like PLRV, many poleroviruses contain a serine protease. Despite their diverse amino acid sequences, these proteases are predicted to share similar structures. Rladg recognizes all ten tested polerovirus proteases, suggesting a conserved structural recognition mechanism. We propose that Rladgs broad recognition capacity could enable resistance to poleroviruses in many crop species. Rladg is the first R-gene reported to confer resistance to a phloem-limited pathogen and could provide enhanced resistance to many economically important poleroviruses.

Rust fungi are significant threats to global food security, causing substantial damage to crops through their ability to adapt and evolve new strains that overcome resistance. These obligate biotrophs infect host plants by secreting effector proteins that manipulate host physiology to promote infection and colonisation. We used AlphaFold2 to investigate structural conservation among effector proteins for the secretomes of Melampsora lini and four Puccinia species. AlphaFold2 yielded high-confidence predictions for 45.7% of the 27,090 secreted proteins, while 19% were poor quality. Comparative analysis revealed extensive structural diversity across the rust secretomes, with all thirteen known rust Avr proteins belonging to different clusters apart from AvrSr13 and AvrSr33. Nevertheless, there were still numerous large clusters of structurally-related proteins, including 59 clusters with over 50 members each, three of which contained known Avr proteins. Of the major structural families defined in other fungi, the rust species studied here only contained FOLD and ToxA-like families. Structural analysis of cysteine-rich proteins revealed over a thousand effector candidates featuring zinc-binding sites, with approximately 75% predicted to be cytoplasmic effectors. In contrast, cysteine-rich apoplastic effector candidates were characterized by a high frequency of disulfide bonds. One family of predicted metal-binding proteins was greatly expanded in P. graminis f. sp. tritici and includes AvrSr13 and AvrSr33. We confirmed that purified AvrSr13 and AvrSr22 proteins bind to zinc in vitro using biochemical assays. Taken together, structural modeling provides new avenues to study sequence-unrelated effectors and highlights the high degree of diversity in the effector repertoires of rust species.

The receptor kinase FERONIA (FER) perceives endogenous RAPID ALKALINIZATION FACTOR (RALF) peptides and regulates a plethora of plant physiological processes, including immunity. RALF peptides also bind to LEUCINE-RICH REPEAT EXTENSIN (LRX) proteins as structural components of the cell wall. We recently showed that the FER-RALF-LRX pathway supports colonization by the obligate biotrophic fungal pathogen Erysiphe cruciferarum (Ecr), a member of the powdery mildew species complex that infects Arabidopsis. Genetic disruption of the pathway primarily affects conidiation of the fungus, raising the question of effects on fungal nutrition. To get further insight into the underlying mechanisms, we performed RNA sequencing (RNAseq) to identify differential transcriptional responses of FER-RALF-LRX pathway mutants upon Ecr infection. Surprisingly, our results revealed that pathway disruption has a limited impact on the overall transcriptional changes upon fungal infection. However, consistent with previous reports, FER-RALF-LRX pathway mutants show changes in basal expression of a plethora of genes, mainly associated with cell wall metabolism, jasmonic acid signalling, amino acid biosynthesis and secondary metabolism. Many of these genes are regulated by Ecr infection across genotypes, too. This raises the question whether these are relevant pathway components for powdery mildew host establishment downstream of the FER-RALF-LRX module. In summary, our data reveals new insights into FER-RALF-LRX-dependent responses that may support host susceptibility to biotrophic plant pathogens.

Trichoderma species are widely used as biological control agents due to their ability to parasitize plant pathogens. However, substantial variability in mycoparasitic performance exists among strains, even within the same species, and the underlying molecular mechanisms remain poorly understood. Here, we performed comparative genomic and transcriptomic analyses of six Trichoderma atroviride strains exhibiting contrasting mycoparasitic performance (weakly or highly parasitic; WP or HP) against Alternaria brassicicola, Rhizoctonia solani, and Globisporangium ultimum. Comparative genomics revealed limited strain-specific differences, mainly restricted to NLR (NOD-like receptor) repertoires, with certain NLR-coding genes absent from WP strain genomes compared to HP strains, while overall genomic variation remained low. In contrast, transcriptomic analyses revealed strong differences in gene expression dynamics between HP and WP strains. Co-expression network analysis identified two modules associated with mycoparasitic performance. The first was specifically induced in response to pathogen contact and was enriched in genes encoding cell wall-degrading enzymes, with stronger expression in HP strains. The second module was more broadly overexpressed in HP strains across all conditions and included genes involved in detoxification and defense-related pathways. In addition, this module encompassed genes involved in specialized metabolite biosynthesis and effector-like protein secretion, with WP and HP strains differentially expressing distinct gene subsets within these categories. Together, these results provide a comprehensive framework for identifying the molecular drivers of mycoparasitic performance in T. atroviride. This study deepens our understanding of the functional diversity within the species and establishes a robust foundation for the future development of molecular markers to predict strain efficiency.

Common scab disease on potato is caused by members of more than 10 pathogenic Streptomyces species. Genome-enabled methods are being increasingly deployed to characterize Streptomyces that cause common scab disease of potato and other tuber and root crops. However, the study of phytopathogenic Streptomyces is constrained by the limited availability of high-quality genome sequences. Here we report improvements to the quality and completeness of genome assemblies for 12 pathogenic type strains of Streptomyces and six closely related non-pathogenic type strains. These assemblies have an average N50 of 7.4 Mbp and with BUSCO scores all greater than 98.5%. Analyses showed that the genomes of phytopathogenic Streptomyces are consistently among the largest Streptomyces genomes sequenced and, relative to those of non-pathogenic strains, are more enriched in genes involved in carbohydrate and amino acid metabolism. Plasmids were not consistently detected across assemblies, suggesting that they are not conserved across species and are not necessary for pathogenicity. Furthermore, comparisons of genome assemblies among both closely and distantly related strains revealed multiple rearrangements within linear chromosomes and reduced synteny near telomeric regions. These improved genome assemblies, many of which correspond to type strains, provide valuable resources for advancing our understanding of the pathogenicity in the genus.

Parastagonospora nodorum (Berk.) Quaedvlieg, Verkley & Crousis, a necrotrophic fungal pathogen, is the causal agent for septoria nodorum blotch, a major constraint on global wheat production. Pathogen-produced necrotrophic effectors (NEs) that interact with host-sensitivity genes in an inverse gene-for-gene manner, collectively leading to effector-triggered susceptibility (ETS). Here, we investigated the transcriptional responses of two Triticum aestivum L. genotypes, Mace and Lancer, following infiltration with a novel NE, SnTox8. A total of 12,679 unique differentially expressed genes in Mace and 149 in Lancer were detected from transcriptomic analysis. In the SnTox8-sensitive cultivar, Mace, numerous defence-related genes were induced, including protein phosphorylation cascades, reactive oxygen species bursts, calcium signalling, phytohormone modulation, and suppression of photosynthesis, consistent with findings from other ETS models, in which necrotrophic fungal pathogens hijack host defence systems to proliferate. The interaction also activated genes involved in signal transduction, metabolism, membrane modification, and molecular transport, reflecting a coordinated host reprogramming that promotes cellular dysfunction and cell death, thereby facilitating necrotrophic pathogenesis. In contrast, Lancer, an SnTox8-insensitive cultivar, exhibited minimal transcriptional changes with no evidence of effector recognition or downstream defence-related activities. Overall, this study exhibited that SnTox8 manipulates kinase-mediated immune signalling and metabolic reprogramming to convert defence activation into host cell death, revealing a mechanistic basis for ETS in wheat. The identified SnTox8-Snn8-triggered processes were confirmed through additional transcriptome analysis of Mace mutants. Outcomes from this study establish a foundation for identifying, functionally characterising and validating the corresponding host susceptibility gene Snn8.

Delivery of therapeutic biomolecules into plant vascular tissues remains a challenge in management of vector-borne plant pathogens. The symbiont concept uses reprogrammed Agrobacterium tumefaciens galls (called symbionts) to produce biomolecules while remaining connected to host vasculature. We evaluated whether symbionts expressing antimicrobial peptides (AMPs) suppress Candidatus Liberibacter solanacearum (CLso), the causal agent of potato zebra chip disease. Symbionts were engineered to express a Streptococcus mutans bacteriocin associated with bacterial membrane disruption (Blp-Sm), or an AMP isolated from finger lime and associated with resistance to citrus greening disease (MaSAMP). Effects of AMP-producing symbionts on CLso titers, infection incidence, pathogen movement, and disease symptoms were evaluated in tomato and potato. In tomato, neither AMP significantly reduced CLso titers or infection incidence. However, in potato, AMP-producing symbionts reduced CLso accumulation and movement from CLso-inoculated source shoots into non-inoculated sink shoots connected through underground tubers. Blp-Sm produced the strongest reduction in CLso accumulation and infection incidence in sink tissues. In separate assays where symbionts were established directly on potato seed tubers, MaSAMP significantly reduced CLso titers in stems and tubers and reduced zebra chip symptoms in tubers, despite no reduction of CLso titers in terminal leaves. These findings demonstrate that AMP-producing symbionts suppress vascular pathogen accumulation and movement within plants and highlight the symbiont concept as a potential platform for managing diseases caused by vascular-restricted pathogens. Further, they show the potato-CLso system is a promising infection model to both refine and improve symbiont technology, and to test additional AMPs for potency against related pathogens.

Many Ceratocystis species cause cankers and unique vascular wilt diseases, often on a broad and unpredictable range of plant hosts. Characteristic necrosis of xylem parenchyma cells and dark staining of surrounding tissue is typically evident, especially in woody hosts. The molecular basis for this unique pathogenicity and host range remains unclear, but bacterial-type phosphatidylinositol phospholipase C (bPI-PLC) genes were recently identified in unusually high copy number in multiple Ceratocystis species, and the PLCs may play a role in host membrane disruption. We produced a high-quality long-read genome assembly of the rapid [O]hia death pathogen, Ceratocystis lukuohia, and identified 81 partial or complete PLC-like genes, each with a unique DNA sequence, encoding signal peptides and a PLC-X domain. The putative translations mostly ranged from 300 to 500 amino acids that differed markedly from the fungal and prokaryotic bPI-PLCs at sites conferring phosphatidylinositol specificity, suggesting a novel family of secreted PLCs (Cer-PLCs). Remarkably, 73 of the full or partial Cer-PLC genes reside in a single 543 kb gene cluster in C. lukuohia. Comparison to an available long-read genome assembly of C. fimbriata revealed a similar Cer-PLC cluster of 61 genes, with a gene order and arrangement broadly similar to that of the C. lukuohia cluster, except for a large inversion at the beginning of the cluster. Differences suggest that the cluster is dynamic, with many apparent indels involving multiple Cer-PLCs. We compared 40 newly-assembled genomes of Ceratocystis strains and eight publicly available genomes and found that the Cer-PLCs comprise a gene family present in all Ceratocystis species but differing greatly in number (26 to 92), with 64 to 92 in species of the highly aggressive Latin American Clade. The two closest relatives of Ceratocystis have Cer-PLCs but not in the gene cluster: Chalaropsis spp. have only one Cer-PLC, and Berkeleyomyces basicola has 25 related Cer-PLCs scattered across multiple contigs. No Cer-PLC was detected in the more-distant members of the Ceratocystidaceae. The unique cluster in Ceratocystis apparently arose through insertion of Cer-PLCs within an ancestral gene cluster with a CeGAL transcription factor, followed by repeated duplications and rapid diversification of Cer-PLCs, perhaps driven by unequal crossover events. This extraordinary expansion, diversification, and maintenance of Cer-PLCs may have played a major role in the evolution of aggressiveness and host range in Ceratocystis. Impact StatementNew strains of Ceratocystis species with expanding host ranges are emerging as important plant pathogens around the world. However, little is known about the basis for the wide variation in host range and aggressiveness of Ceratocystis species. An earlier study had identified a gene family coding for phosphatidylinositol-specific phospholipase C (PI-PLC) in some Ceratocystis species. Our sequence analyses suggest that the coded enzyme is not likely phosphatidylinositol-specific but may have retained capability of degrading plant membranes and may be a major determinant of aggressiveness and host range. The most aggressive species in the genus has up to 92 copies of this unique class of PLCs, defined here as Cer-PLCs, making the expansion of this gene family among the largest known in fungi. Most of the Cer-PLC genes occur in a gene cluster of more than 500 kb, which appears to be under the control of a CeGAL-type transcription factor. Coordinated regulatory control may enable the hyper-production of these membrane-degrading enzymes during pathogenesis. The Cer-PLC gene family occurs in close relatives of Ceratocystis, but the Cer-PLC gene cluster is unique and universal in Ceratocystis. The gene cluster is the largest known for a single gene family, and it is highly dynamic and likely undergoes frequent recombination. Multiple introductions of strains to a new environment could generate very aggressive recombinants that attack previously unrecognized hosts, as appears to be happening with the multiple introductions of the South American species C. manginecans to Asia.

Blumeria graminis f. sp. avenae (Bga), the causal agent of oat powdery mildew, is one of the most host-specialized members of the B. graminis species complex. Despite its agricultural importance, the lack of a high-quality reference genome has limited studies of host specialization, virulence evolution and comparative genomics in this pathogen. Here, we generated the first chromosome-scale genome assembly of Bga using an integrative approach combining long- and short-read sequencing, Hi-C scaffolding and transcriptome data. The Bga genome exhibits hallmark features of powdery mildew fungi, including extensive repeat content and low gene density. Comparative analyses revealed that genome expansion is primarily associated with historical transposable element proliferation rather than recent transpositional activity. Genome organization is consistent with a functionally stratified "one-speed" model, in which genes associated with pathogenicity, including predicted effectors and infection-responsive genes, are preferentially located in transposable element-rich regions characterized by reduced synteny conservation and extended intergenic spaces. In contrast, conserved genes are concentrated in compact genomic regions and maintain strong syntenic conservation across cereal-infecting formae speciales. Hi-C analyses demonstrated a highly structured chromatin architecture and revealed genome organization patterns associated with infection-related gene expression. Comparative genomic analyses indicated that host specialization in Bga is driven by localized diversification of a relatively small subset of genes rather than large-scale genome restructuring. These results provide the first high-quality genomic resource for Bga and offer new insights into the evolutionary mechanisms underlying host specialization in powdery mildew fungi.

Symbiotic interactions between plants and nitrogen-fixing microorganisms are essential for sustainable agriculture, yet the molecular mechanisms underlying plant-cyanobacterium symbiosis remain poorly understood. In particular, the nature of the signalling mechanisms mediating partner recognition in associations involving Nostoc species is largely unknown. Recent proteomic analyses have identified proteins homologous to rhizobial Nod factors biosynthetic enzymes in Nostoc punctiforme, suggesting the existence of a Nod-like signalling system. However, the functional role of these components has not been experimentally validated. Here, we investigate the contribution of nod-like biosynthetic and regulatory genes to symbiosis by analysing mutants of N. punctiforme affected in genes with homology to nodB and nodD. Phenotypic characterization revealed that disruption of nodB-like genes does not impair free-living growth but affects early stages of plant association and colonization. Specifically, the nodB1 mutant is impaired in plant association and shows a mild defect in colonization, whereas the nodB3 mutant exhibits a severe defect in colonization. In contrast, nodD-like mutants exhibited altered symbiotic phenotypes, with specific regulators differentially affecting interaction and colonization efficiency in rice (Oryza sativa). In particular, mutation of nodD2 and nodD3 reduced plant association and severely compromised colonization in Oryza sativa, with a more pronounced phenotype in nodD3 mutant. Altogether, our results provide genetic evidence supporting the involvement of Nod-like components in cyanobacterial symbiosis and suggest the existence of a regulatory and biosynthetic module contributing to plant colonization. These findings shed new light on the evolution and diversity of symbiotic signalling mechanisms across plant-microbe interactions.

Tomato worldwide production is increasingly challenged by complex disease outbreaks involving multiple interacting pathogens. In Israel, recent years have seen a marked rise in vascular collapse symptoms in greenhouse-grown tomatoes, coinciding with the widespread emergence of tomato brown rugose fruit virus (ToBRFV) and pepino mosaic virus (PepMV). Here, we investigated the bacterial and viral agents associated with these outbreaks and examined how viral infection influences the development and severity of bacterial diseases that cause vascular collapse. Surveys conducted between 2021 and 2026 revealed that tomato pith necrosis outbreaks were associated with a diverse bacterial community dominated by members of the Pseudomonadales and Enterobacterales, while bacterial canker outbreaks were exclusively linked to Clavibacter michiganensis. Multilocus sequence analysis showed that pith necrosis-associated Pseudomonas isolates clustered primarily within the P. syringae, and P. corrugata phylogroups. Pathogenicity assays demonstrated that only a subset of pith necrosis-associated bacteria, P. mediterranea, P. capsici, P. viridiflava, and Xanthomonas euvesicatoria pv. perforans, induced pith necrosis under controlled conditions, with high variability in symptom severity. Co-inoculation experiments showed that ToBRFV- and PepMV-infected plants exhibited a 40-100% increase in lesion size following inoculation with pith necrosis-associated bacteria, without a corresponding increase in bacterial colonization, whereas the same viral infections attenuated wilt symptoms caused by C. michiganensis. Together, our findings demonstrate that endemic viral infections differentially modulate bacterial disease outcomes, either exacerbating or attenuating symptoms depending on the pathogen. These results highlight the importance of multi-pathogen interactions in disease severity and have important implications for tomato disease management.

Ada1 (All Development Altered-1) is a conserved but poorly characterised basic leucine zipper (bZIP) transcription factor found throughout filamentous fungi. In the wheat pathogen Parastagonospora nodorum, PnAda1 is required for full virulence and is transcriptionally associated with the virulence regulator PnPf2, but its biological functions remain unclear. Here, we combined comparative RNA sequencing with targeted phenotypic analyses to define the role of PnAda1 during vegetative growth and host infection. Deletion of PnAda1 did not abolish pathogenicity but delayed disease progression, with the PnAda1-deletion mutant transcriptome at 7 days post-inoculation resembling that of the wildtype SN15 at 3 days. This developmental delay was associated with impaired activation of early infection-associated genes, including putative carbohydrate-active enzymes, proteases, transporters and other host-colonisation factors. In contrast, expression of major necrotrophic effector genes was not reduced and instead remained elevated during later stages of infection, indicating that PnAda1 is required for the timely progression of infection-associated transcriptional regulation rather than direct activation of effector genes. Beyond virulence, transcriptomic and phenotypic analyses revealed roles for PnAda1 in nitrogen assimilation, carbon utilisation, abiotic stress responses and fungicide sensitivity. Notably, PnAda1 deletion increased sensitivity to succinate dehydrogenase inhibitor fungicides and reduced expression of succinate dehydrogenase subunit genes. Collectively, our findings identify PnAda1 as a broad regulator of developmental and infection-associated transitions in P. nodorum and expand current understanding of the transcriptional network underlying virulence, metabolism and stress adaptation in an important fungal wheat pathogen.